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Abstract 

In this document, the first three years of a time series of bio-optical marine and atmospheric measurements 
are presented and analyzed. These measurements were performed from an oceanographic tower in the northern 
Adriatic Sea within the framework of the Coastal Atmosphere and Sea Time Series (CoASTS) project, an ocean 
color calibration and validation activity. The data set collected includes spectral measurements of the in-water 
apparent (diffuse attenuation coefficient, reflectance, Q-factor, etc.) and inherent (absorption and scattering 
coefficients) optical properties, as well as the concentrations of the main optical components (pigment and 
suspended matter concentrations). Clear seasonal patterns are exhibited by the marine quantities on which an 
appreciable short-term variability (on the order of a half day to one day) is superimposed. This short-term 
variability is well correlated with the changes in salinity at the surface resulting from the southward transport 
of freshwater coming from the northern rivers. Concentrations of chlorophyll a and total suspended matter span 
more than two orders of magnitude. The biooptical characteristics of the measurement site pertain to both 
Case-1 (about 64%) and Case-2 (about 36%) waters, based on a relationship between the beam attenuation 
coefficient at 660 nm and the chlorophyll a concentration. Empirical algorithms relating in-water remote sensing 
reflectance ratios and optical components or properties of interest (chlorophyll a, total suspended matter, and 
the diffuse attenuation coefficient) are presented. 


1. INTRODUCTION 

The use of remotely-sensed ocean color data has in- 
creased considerably over the last years, since several new 
sensors were put in orbit more than 10 years after the 
Coastal Zone Color Scanner (CZCS) ceased its activity. 
These sensors, specifically, the Sea-viewing Wide Field-of- 
view Sensor (SeaWiFS), the Ocean Color and Tempera- 
ture Scanner (OCTS), the Marine Optical Spectroradiome- 
ter (MOS), the Moderate Resolution Imaging Spectrora- 
diometer (MODIS), and the Medium Resolution Imaging 
Spectrometer (MERIS), present advanced characteristics 
compared to CZCS, in the addition of a certain number 
of channels, and an increased signal- to- noise ratio and dy- 
namic range (IOCCG 1998). These increased capacities 
are intended, in particular, to improve the retrieval of chlo- 
rophyll a concentration in the open ocean (Hooker et al. 
1992 and O’Reilly et al. 1998), as well as eventually allow- 
ing the retrieval of additional marine quantities such as 
total suspended matter or colored dissolved organic mat- 
ter (Carder et al. 1999). 

For the open ocean, the possibility of retrieving chlo- 
rophyll a and additional water components from remotely- 
sensed data is presumably simplified, because in most cases, 
they all co-vary (but see Carder et al. 1991). In coastal wa- 
ters (where most of the so-called optical “Case-2” waters 
can be found), however, this retrieval is less successful. 
In effect, because of the diversity of the sources (rivers, 
atmospheric input, bottom resuspension, etc.) the sim- 
plifying assumption of co-variation may not be verified. 
Consequently, by contributing — often antagonistically — to 
the total marine signal, the occurrence of these different 
components makes their separation much more difficult on 
an optical basis. 


As a consequence of the number of sensors now opera- 
tional, there is an increasing demand for in situ bio-optical 
data (water and atmosphere) in order to address two main 
topics: 

1) The development of models and algorithms allowing 
the retrieval of marine and atmospheric quantities 
from the signal at the top of the atmosphere; and 

2) The vicarious calibration of the sensors, and the 
validation of their different products. 

Programs including time-series of measurements on spe- 
cific sites such as the Marine Optical Buoy (MOBY, Clark 
et al. 1997), the Yamato Bank Optical Mooring (YBOM, 
Kishino et al. 1997), and the Plymouth Marine Bio-Optical 
Data Buoy (PlyMBODy, Pinkerton and Aiken 1999), or 
repetitive ship cruises such as the Atlantic Meridional Tran- 
sect (AMT) program (Aiken et al. 2000) were developed 
considerably in the past years. Basing the investigation 
of the environment on a long time or space scale, one of 
their objectives was, in particular, to encompass a wide 
range of optical properties and concentration of optically 
significant components. 

In this context, a program of extensive in situ mea- 
surements called the Coastal Atmosphere and Sea Time 
Series (CoASTS), was set up (Zibordi et al. 2002). Within 
this program, a series of marine and atmospheric measure- 
ments are performed from an oceanographic tower located 
in the northern Adriatic Sea, off the Venice Lagoon. As of 
this writing, the project has provided more than six years 
of data (with a frequency of 1-2 campaigns per month) 
including water apparent and inherent optical properties 
and the corresponding optically significant components, as 
well as atmospheric measurements, which allows for the es- 
timation of aerosol properties. The resulting data set was 
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partially presented (Berthon et al. 1998, 2000, 2001; and 
Zibordi and Berthon 2001) and used for radiative transfer 
computations (Zibordi et al. 1999), and calibration and 
validation activities (Sturm and Zibordi 2002). 

This report is aimed at presenting the CoASTS data 
from October 1995 to December 1998, in particular by giv- 
ing emphasis to a number of items: 

1 ) The variability at different time scales of some rel- 
evant quantities (the main optical components and 
their inherent optical properties, the apparent opti- 
cal properties of the water column, etc.); 

2 ) Inter-relationships between quantities to be used in 
the development of bio-optical models; and 

3) Local algorithms relating the concentration of to- 
tal chlorophyll a (C a ) and the concentration of total 
suspended matter (Ctsm)- 

2. DATA AND METHODS 

Within the framework of the CoASTS project, the Joint 
Research Centre with the contribution of the Italian Na- 
tional Research Council, has sponsored measurements that 
are performed from the Acqua Alta Oceanographic Tower 
(AAOT), located in the northern Adriatic Sea, off the 
Venice Lagoon (45.310°N,12.510°E). Measurement cam- 
paigns lasting a few days, take place every 2-4 weeks ac- 
cording to the season and can include from 3-10 sets of 
comprehensive measurements (or stations) each. 

The measurements conducted during the periodic field 
campaigns include in-water profiles of the apparent and 
inherent optical properties of the light field, discrete water 
samples at different depths (near-surface, 8 m, and 14 m; 
the bottom is at 17 m) for biogeochemical analysis, temper- 
ature and conductivity profiles, sky radiance, sun irradi- 
ance, and meteorological data. The CoASTS data set pre- 
sented here includes 269 stations (October 1995-December 
1998). A detailed description of the instruments, measure- 
ment methodologies, and uncertainties is given in Zibordi 
et al. (2002). 

3. DATA ANALYSIS 

In this section, the time series of the marine and at- 
mospheric measurements are presented and discussed. 

3.1 Environmental Characteristics 

The surface circulation in the northern Adriatic Sea 
is dominated by the North Adriatic (NAd) current which 
flows southward along the Italian northeastern coast, with 
a maximal southward extension during the boreal spring. 
In the boreal autumn, the surface circulation is primarily 
determined by the North Adriatic cyclonic gyre, which de- 
velops at the latitude of the Po River outflow (Artegiani 
et al. 1997b). Hydrological features at the tower site are 


primarily influenced by the northern river discharges (con- 
veyed by the NAd current), as well as by wind and rain 
variability, and consequently are susceptible to very short 
time-scale variations, which are superimposed on the main 
seasonal evolution. 

3.1.1 Meteorological Conditions 

In Table 1, data are presented of the average seasonal 
values of some meteorological variables (air temperature, 
relative humidity, atmospheric pressure, wind speed, and 
direction) instantaneously measured at the platform site at 
the time of sampling. The use of the calendar months for 
partitioning the seasons is in agreement with Zavatarelli 
et al. (1998) and it maintains coherence with the usual 
seasonal cycle of river discharge. A direct comparison can, 
therefore, be made with the climatology of biogeochemical 
variables presented in the latter work. 

Air temperature, T a , exhibits a total variation of 25° C 
with a minimum around 3°C in January 1998 (an excep- 
tional unique value of 0.3° C was observed in January 1996) 
and a maximum around 28° C in August 1998. The ob- 
served seasonal cycle is in good agreement with the mean 
climatology presented by Artegiani et al. (1997a). Rel- 
ative humidity (RH)> although highly variable at short 
time scale, presents average maximum values in springf 
(April-May-June) of 78.8%, and minimum values in win- 
ter (January-February-March) of 61.7%. 

Wind speed (W 3 ) and wind direction (' Wd ) measured at 
the time of sampling are highly variable and no clear cy- 
cle is discernible. The strongest and lightest winds, how- 
ever, are generally observed during the winter and autumn 
(October-November-December), respectively. Northeast- 
erly (0-90°) and southeasterly (90-180°) winds generally 
dominate, but abrupt changes in wind direction have fre- 
quently been observed during a field campaign. No partic- 
ular correlation has been observed between wind speed and 
direction. Because measurements are not performed when 
weather conditions are too severe, which is usually associ- 
ated with low atmospheric pressure (P a ), the time series 
of meteorological data may not be fully representative of 
the site conditions during the year. 

3.1.2 Near-Surface Temperature and Salinity 

The near-surface values of seawater salinity S w and 
temperature T w presented here correspond to the average 
value within the near-surface layer (for a depth ranging 
from 1-2 m) and computed from conductivity, tempera- 
ture, and depth (CTD), vertical profiles performed during 
sampling stations. Density excess, 7, was computed ac- 
cording to the UNESCO (1981) formulation without in- 
cluding pressure effects. Average global and seasonal val- 
ues for these quantities are presented in Table 1. 


f References throughout this document to seasons of the year 
are to the seasons of the Northern Hemisphere and are de- 
fined on first usage. 
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Tkble 1. Marine and atmospheric parameters for the entire (global) data set and a seasonal partitioning of the 
data: January-February-March (JFM), April-May-June (AMJ), July- August-September (JAS), and October- 
November-December (OND). For each variable (X), X n is the number of samples, X is the average, X a is the 
standard deviation, X~ is the minimum value, and X+ is the maximum value. 


Period 

Statistic 

Ta [°C] 

RH [%} 

P a [hPa] 

W s [ms x ] 

w&[°] 

T w [°C] 

[PSU] 

7 [kgm 3 ] 

Global 

a„ 

259 

265 

268 

247 

238 

269 

269 

269 


A 

15.84 

71.1 

1020.16 

3.31 

144.0 

16.907 

34.978 

25.373 


x„ 

6.82 

17.0 

6.87 

1.90 

107.5 

5.913 

1.774 

2.326 


X- 

0.3 

24 

1000.5 

0.0 

0 

7.31 

30.54 

19.67 


A + 

28.2 

100 

1036.7 

9.1 

360 

27.36 

38.05 

29.60 

JFM 

x n 

56 

63 

63 

63 

62 

64 

64 

64 


X 

7.82 

61.7 

1024.31 

3.74 

137.7 

9.740 

36.121 

27.870 


x« 

2.89 

19.4 

8.19 

1.89 

121.5 

0.914 

1.836 

1.447 


X- 

0.3 

24 

1002.2 

0.0 

0 

7.31 

31.93 

24.43 


a + 

13.9 

99 

1036.7 

8.6 

360 

11.66 

38.05 

29.60 

AMJ 

A n 

68 

65 

68 

63 

61 

68 

68 

68 


X 

17.39 

78.8 

1018.76 

3.26 

144.1 

17.816 

34.375 

24.759 


At 

4.16 

11.6 

6.87 

2.02 

84.3 

4.402 

1.475 

1.756 


x~ 

8.4 

57 

1000.5 

0.0 

5 

11.18 

30.89 

20.33 


A + 

26.0 

99 

1035.4 

8.4 

332 

24.84 

37.28 

28.52 

JAS 

X n 

61 

63 

63 

47 

43 

63 

63 

63 


X 

23.63 

74.9 

1018.19 

3.55 

119.3 

24.283 

34.258 

22.988 


x„ 

2.38 

17.2 

4.21 

2.14 

86.9 

1.731 

1.574 

1.545 


X- 

17.5 

24 

1004.7 

0.0 

6 

20.28 

30.54 

19.67 


A + 

28.2 

99 

1025.0 

9.1 

341 

27.36 

36.89 

25.44 

OND 

An 

74 

74 

74 

74 

72 

74 

74 

74 


A 

14.06 

68.9 

1019.59 

2.82 

164.1 

15.991 

35.156 

25.809 


x a 

5.67 

14.6 

6.08 

1.51 

120.8 

3.702 

1.607 

1.513 


x~ 

3.1 

37 

1008.0 

0.0 

0 

8.20 

30.67 

23.12 


X* 

23.0 

100 

1031.2 

8.0 

360 

21.01 

37.60 

28.66 


The seasonal cycle for near-surface temperature (Fig. 1) 
closely follows the air temperature, with the minimum 
(about 6-8° C) in January-February and the maximum 
(around 26-28° C) in July-August. This cycle is in good 
agreement with the seasonal climatology provided by Arte- 
giani et al. (1997b). Minor intra-campaign changes (around 
1°C) are usually observed except for campaigns 21 (10- 
7.5°C), 37 (11-9°C) and 44 (22.5-24.5°C) for which vari- 
ations larger than 2°C were recorded. These short-term 
variations are frequently observed in conjunction with the 
presence of less saline water transported to the site. The 
latter is usually seen as a thin (about 2m thick), generally 
cooler, surface layer. A variation of about 2°C (from 21- 
19°C), due to the mixing of the water column induced by 
high winds, can also be observed for campaign 11. In this 
last case, near-surface salinity increases from 33.5-35 PSU 
after mixing with more saline subsurface waters. 

High salinity at the surface (Fig. 1), with a maximum 
around 38 PSU, is systematically associated with the win- 
ter months, whereas minimum values are irregularly dis- 
tributed within the year. Very low values, from 30-31 PSU, 
were observed in October 1996 (campaign 20) and, to a 


lesser extent, in July 1997 (campaign 31) and June 1998 
(campaign 44). Note that anomalous high air temperature 
(and low atmospheric pressure) were also observed during 
campaign 20. 

The 1996 campaign year departed from the mean cli- 
matology proposed by Artegiani et al. (1997b), by globally 
showing higher values (differences up to 4-5 PSU), whereas 
years 1997 and 1998 more closely followed it. Rapid chang- 
es, on the order of 2-4 PSU in the near-surface salinity, 
were frequently observed.! These changes axe primarily 
due to variations in the extension toward the east of the 
southward flowing coastal current carrying low salinity wa- 
ters from the northern rivers. 

3.1.3 Vertical Distribution of T w and S w 

The average seasonal vertical profiles of temperature 
and salinity are presented in Fig. 2. In winter (Figs. 2a 
and b), cold waters extend from the surface down to the 


t These changes were seen during campaigns 8, 11, 13, 21, 22, 
30, 31, 36, 37, 39, 41, 43, 44, and 50. 
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Fig. 1 . Time-series of near-surface a) temperature, T w ; and b) salinity, S w . For each campaign, the average, 
the standard deviation around the average, and the individual values are represented by solid circles, a black 
line, and open circles, respectively. The numbers above the symbols refer to the campaign number. 


bottom (the density excess 7 can reach values up to 29.4- 
29.6 characteristic of the northern Adriatic Sea deep wa- 
ters. Nevertheless, the average salinity profile shows lower 
values at the surface (a difference of about 1 PSU with re- 
spect to the values at 5 m). This inhomogeneity is due 
to short time scale events showing the presence of a thin 
(about 2-3 m thick) neax-surface layer of low salinity (gen- 
erally cooler) transported from the north, e.g., campaigns 
6 (January 1996), 21 (January 1997), 22 (February 1997), 
37 (December 1997), or 39 (February 1998). This fact is 
also highlighted by the elevated variability (standard de- 
viation) around the average surface value (Fig. 2b). 

In the spring (Figs. 2c-d), the increase in the air tem- 
perature and the increase in river discharges leads to a 
more stable stratification (e.g., campaign 8 in late March 
1996, or campaign 23 in early April 1997). The thickness 
of the (homogeneous) near-surface layer remains variable, 
probably according to the magnitude and duration of wind 
events. Temperature, in general, progressively decreases 
from the surface to depth, whereas salinity shows an in- 
verse profile. In summer (Figs. 2e-f), a warm (22.5-26°C 
at 1 a) and low salinity (32.5-36.0 PSU) layer extends from 
the surface down to about 10 m. At the bottom, tempera- 
ture and salinity reach 18-22° C and 36.5-37.5 PSU, respec- 
tively. Here also, salinity and temperature show markedly 
opposite vertical profiles. The autumn season (Figs. 2g-h) 
displays the opposite influence of autumnal river discharges 
and cooling of the air temperature: salinity and tempera- 
ture both increase, more or less abruptly, from the surface 
to depth. 


This average or generalized description of the seasonal 
evolution is in good agreement with the one presented for 
the whole “Shallow Northern Adriatic” Sea by Zavatarelli 
et al. (1998). 

3.2 Atmospheric Optical Characteristics 

CoASTS atmospheric measurements include measure- 
ments of the sky radiance in the sun and the almucantar 
planes, and direct sun irradiance, as well as the ratio of 
diffuse-to-direct solar irradiance, both continuously mea- 
sured. The data presented in this section axe the aerosol 
optical thickness, r a (A), where A is the wavelength, and a 
is the Angstrom exponent. The wavelength dependency of 
the aerosol optical thickness is commonly expressed by the 
Angstrom law (Angstrom 1961) as: 

r a (A) = /?A"“, (1) 

where 0 is the Angstrom coefficient. The r a values were 
obtained from measurements performed on the tower by a 
sun photometer included in the Aerosol Robotic Network 
(AERONET) network (Holben et al. 1998) managed by the 
National Aeronautics and Space Administration (NASA) 
Goddard Space Flight Center (GSFC). The r a data were 
screened for clouds by removing any data exceeding the 
threshold value of 0.1 for r a (870) and then only considering 
situations for which the relative air mass m was less than 
or equal to 5. Values of a and 0 were then computed by 
least-squares fit of the estimated r a ( A) versus A for the 
wavelengths 440, 670, and 870 nm. 
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Fig. 2. Average seasonal vertical profiles (lm binning) of temperature [panels a) JFM, c) AMJ, e) JAS, g) 
OND] and salinity [panels b) JFM, d) AMJ, f) JAS, h) OND]. Standard deviations around the average are 
represented by dotted lines; n refers to the number of stations used for computing the averages. 
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Fig. 3. Time-series of a) the aerosol optical thickness r 0 , at 555 nm, and b) the Angstrom exponent a. Filled 
circle and black lines represent the daily averages and standard deviations, respectively. 


The aerosol optical thickness r a at 555 nm (computed 
using the retrieved values of a and (3) and the exponent a, 
are both shown in Fig. 3 as daily average values together 
with the associated standard deviation. Values of r a (555) 
vary in the range of 0.02-0.26 with a global average value 
of 0.136 (a = 0.055). 

Standard deviations of daily average values (Fig. 3a), 
indicate significant short-term variations of the aerosol op- 
tical thickness. These can be attributed to short-term 
changes in the aerosol composition and concentration, due 
to changes in the local circulation. In Fig. 3b, the exponent 
a, varies from 0.7-2. 4, and shows an average value of 1.70 
(a = 0.27). These values of a are in agreement with inde- 
pendent measurements performed at different sites in the 
Po River valley (De Santis et al. 1994) and indicate that 
small aerosol particles prevail over large particles, thus sug- 
gesting a general prevalence of continental aerosol in the 
measurement area. 

3.3 Marine Optical Components 

Marine optical components represent the different el- 
ements within the water column (including water itself) 
contributing to the attenuation of light either by absorb- 
ing or by scattering photons. Three significant compo- 
nents are classically considered for bio-optical modeling 
in coastal areas: the pigmented particulate matter (phy- 
toplankton pigments), the nonpigmented particulate mat- 
ter, and the colored dissolved organic matter — also called 


“yellow substance” — expressed through its absorption co- 
efficient (presented later in the section on inherent optical 
properties, Sect. 3.4). 

3.3.1 Phytoplankton Pigments 

The concentration of phytoplankton pigments (Table 2) 
was measured by high performance liquid chromatogra- 
phy (HPLC, JGOFS 1994). The total chlorophyll a (C a ) 
concentration for near-surface samples varies from 0.098- 
8.255 mg m -3 with a global average of 1.309 (a = 1.141, 
Table 2). Depth-integrated concentrations (from 1-14 m) 
vary from 1.7-64.6 mgm -2 . In Fig. 4a, a rough seasonal 
cycle can be observed with a first maximum in February- 
March (eventually May) and a second one in October- 
November; minimum values are observed in December- 
January and June-July. The minima are associated with 
periods of low vertical stability in December- January, or 
low nutrient content-stratified waters in summer; whereas 
the maxima generally correspond to moderately stratified, 
nutrient-enriched, waters of late winter-spring and autumn 
(Gilmartin et al. 1990 and Zavatarelli et al. 1998). 

Considering the average values for the four seasons (Ta- 
ble 2), the cycle is in agreement with the one defined by 
Zavatarelli et al. (1998) for the shallow northern Adri- 
atic Sea (north of the Po River delta): maximal in winter 
and autumn, minimal in summer. Gilmartin and Reve- 
lante (1980) observed three different maxima of C a (fluori- 
metric method) in January-February, May- June- July, and 
October-November; however, this observation was made 
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Table 2. Surface pigment concentrations (in milligrams per cubic meter), including chlorophylls c\ and C2, 
diadinoxanthin, alloxanthin, diatoxanthin, and /3-carotenoids — C\ 2, C^, Ca, Cd, and Cp, respectively— and 
total suspended matter concentrations, Ctsm in grams per cubic meter. The period and statistic notations are 


the same as in 


Table 1. 


Period 

Stat. 

Ca 

C b 

C 12 

Cp 

Chb 

Cp 

C A 

Cd 


Cfi 

Ctsm 

Global 

x n 

268 

259 

233 

267 

268 

265 

259 

243 

263 

249 

268 


X 


0.0941 

0.3005 

0.4293 


EBSI 



0.0877 

0.0531 

1.081 


Xo 


0.1155 

0.4468 

0.5520 

0.1791 

H 

Hi 

0.0576 

0.0944 

0.0538 

0.695 


X- 


0.007 

0.008 

0.008 


Hi 

Hi 

mm 

0.004 

0.002 

mm 


X* 

8.255 

0.763 

3.022 

3.267 

1.163 

EH 

EH 

EH 

0.608 

0.363 

H 

JFM 

X n 

64 

58 

51 

64 

64 

61 

60 

48 

59 

59 

64 


X 

1.5930 

0.1017 

0.5837 

0.6723 

0.2045 

0.1384 

0.0584 

0.0724 

0.0283 

0.0481 

0.992 


X„ 

1.6860 

0.0886 

0.8419 

0.8719 

m 

0.1616 

0.0703 

0.0752 

0.0128 

0.0501 



X- 

0.098 

0.011 

0.008 

0.008 

1 

0.002 

0.002 

0.001 

0.005 

0.003 

0.27 


x + 

8.255 

0.293 

3.022 

3.267 

1.163 

0.670 

0.347 

0.242 

0.062 

0.166 

2.33 

AMJ 

Xn 

68 

68 

68 

68 

68 

68 

67 

65 

68 

59 

68 


X 

1.1494 

0.0485 

0.2118 


0.3309 

0.1306 

0.0516 

0.0645 

0.0736 

0.0459 



■ X„ 

0.9553 

0.0248 

0.1963 

0.3686 

0.1936 

0.1457 

0.0451 

0.0562 

0.0564 

0.0435 



x~ 

0.194 

0.010 

0.016 

0.012 

0.053 

0.007 

0.003 

0.002 

0.004 

0.002 

0.20 


X* 

4.964 

0.101 

0.880 

1.799 

0.925 

0.720 

0.197 

0.236 

0.246 

0.200 

4.13 

JAS 

Xn 

63 

60 

59 

63 

63 


60 

63 

63 

58 

62 


X 

0.9278 

0.0369 

0.1488 

0.3055 

0.1746 


0.0166 

0.0439 

0.1333 

0.0538 

1.200 


■ X, 

0.6172 

0.0279 

0.0976 

0.2789 

0.1091 


0.0121 

0.0565 

0.0861 

0.0396 

0.940 


x~ 

0.186 

0.007 

0.018 

0.035 

0.031 


0.003 

0.005 

0.015 

0.004 

0.27 


X* 

3.740 

0.165 

0.563 

1.778 

0.444 


0.051 

0.346 

0.414 

0.189 

5.07 

OND 

Xn 

73 

73 

55 

72 

73 

73 

72 

67 

73 

73 

74 


X 

1.5391 

0.1774 

0.3106 

0.4629 

0.2477 

0.0973 

0.1018 

0.0265 

0.1095 

0.0624 

0.953 


x a 

0.9406 

0.1691 

0.1796 

0.4258 

0.1636 

0.0852 

0.1116 

0.0300 

0.1327 

0.0708 

0.402 


X~ 

0.336 

0.007 

0.034 

0.051 

0.009 

0.009 

0.006 

0.002 

0.010 

0.003 

0.27 


. x + 

5.163 

0.763 

0.811 

1.807 

0.698 

0.481 

0.447 

0.187 

0.608 

0.363 

2.00 


within a region of the western part of the northern Adri- 
atic Sea which is mainly influenced by the Po River dis- 
charge. Appreciable short time-scale variability adds to 
the seasonal background. 

On several occasions, the value of C a changed within 
a campaign (and often within 24 h) by a factor of 2-3 
(during October 1995, May 1996, January-February 1997, 
October-November 1997, May 1998, and October 1998), 
and up to 8-10 (February-March 1996, and July 1997). 
This short time-scale variability was often observed in con- 
junction with a rapid lowering of the near-surface salinity 
which indicates the southward transport of material from 
the discharge of northern rivers. In particular, an increase 
of C a from 50-850% (computed for pairs of successive sta- 
tions) were observed during situations of decreasing salin- 
ity (a decrease of 0 to — 3PSU). These events can also lead 
to the observation of C a maxima in late winter (February 
or March). Situations of decreasing C a corresponding to a 
decrease in salinity, do not lead to a decrease of C a higher 
than 50%. On the contrary, incidents of increasing C a , to- 
gether with an increase in salinity, were very few and were 
not found to be significant. 


Pigments, particularly the carotenoids, are known to be 
taxonomic markers of phytoplanktonic populations. The 
concentration of chlorophyll b (C*>), which is an indicator 
of Chlorophycea (green flagellates), shows maximum val- 
ues in December or January (Fig. 4b), i.e., when the sea- 
water is cold and well mixed. The latter is in agreement 
with the observation that green flagellates seem to require 
strong mixing situations to develop (Bustillos-Guzman et 
al. 1995). Chlorophyll b concentrations can range from 
0.007-0.763 mgm“ 3 . 

The fucoxanthin concentration (Cp), which is an indi- 
cator of the presence of diatom populations, varies from 
0.008-3.267 mg m -3 . This closely follows the C a cycle, be- 
cause the diatom populations are important contributors 
to the total C a (Fig. 4c). Diatoms are considered to be very 
opportunistic; they optimally develop during transition 
phases between mixed and stratified conditions (Claus- 
tre et al. 1994 and Bustillos-Guzman et al. 1995). Their 
maxima can, therefore, be observed in spring and autumn. 
Surface stability and nutrient resources available in some 
situations (surface transport of lower salinity-nutrient en- 
riched waters) are also favorable to their growth. 
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Fig. 4. Time-series of the concentrations of a) total chlorophyll a (C a ), b) chlorophyll b (C^), c) fucoxanthin 
(CV), and d) the sum of the concentrations of 19'-hexanoyloxyfucoxanthin and 19'-butanoyloxyfucoxanthin 
{Chb)- The symbols and lines are the same as presented in Fig. 1. The values represent the average for each 
campaign. 
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Fig. 4. (cont.) Time-series of the concentrations of e) zeaxanthin (C z ). The symbols and lines are as in 
Fig. 1. Panel f) shows the time-series of the ratio C/C a where C is the pigment, which in this case is C\> 
(dotted line), Cp (solid line), Chb (dashed line), and Cz (dotted and dashed line). The values represent the 
average for each campaign. 


The concentration of 19'-hexanoyloxyfucoxanthin and 
19'-butanoyloxyfucoxanthin combined (C/m), which is an 
indicator of prymnesiophyte and chrysophyte populations, 
ranges from 0.009-1.163 mg m“ 3 and follows the fucoxan- 
thin cycle with a one month shift (Fig. 4d). These popu- 
lations are known to be able to adapt to a variety of water 
column structures (Bustillos-Guzman et al. 1995). 

The zeaxanthin concentration, C z (which is an indi- 
cator of cyanobacteria and prochlorophyte populations), 
ranges from 0.004-0.608 mg m~ 3 . The seasonal variations 
of C z clearly show a shift in time with respect to fuco- 
xanthin concentrations, and are almost perfectly opposite 
to the Cb cycle, with a very distinct maximum in August- 
September when the water is thermally stratified, and a 
minimum in January-February when the water is cold and 
well mixed (Fig. 4e). This global succession can be high- 
lighted by the C:C a ratios (Fig. 4f): 

■ The Ct,:C Q ratio reaches maximum values of 0.2-0. 4 
in winter; 

■ The Cp'C a ratio reaches values of about 0.4-0. 7 in 
March-April or October-November; 

■ The Chb'-Cq, reaches values of 0.35-0.65 mainly in 
April-May or November; and 

■ The C z :C a ratio reaches values of 0.25-0.35 in June- 
July-August. 

This temporal succession is in agreement with previous 
analyses reported in the literature, particularly for the 
Mediterranean Sea (e.g., Bustillos-Guzman et al. 1995). 


Concentrations of the other measured pigments can be 
found in Table 2. 

3.3.2 Total Suspended Matter 

The concentration of total suspended matter, Ctsm> 
was measured according to Strickland and Parson (1972). 
Values ranged from 0.20-5.07 gm’ 3 with a global average 
of 1.081 (a— 0.695), as shown in Table 2. The basic sea- 
sonal cycle of Ctsm (Fig. 5a) shows maxima in April-May 
and October-November and minimum in January. Nev- 
ertheless, the second minimum observed in July-August 
for C a is not as well evidenced here. The short time-scale 
variability characterizing C a also affected the Ctsm dis- 
tribution near the surface: a decrease in salinity of —0.5 
to -2.5 PSU was related to a Ctsm increase of from 100- 
700%. Such events were particularly important during the 
summers of 1997 and 1998, and led to a global average 
summer concentration equivalent to the one in the spring 
but with a higher standard deviation (Table 2). It should 
be noted that the concentrations from summer 1997 (in 
particular, campaign numbers 30 and 31) were affected by 
the presence of high quantities of suspended mucus mate- 
rial as part of a so-called mucilage event, which occurred 
in the northern Adriatic Sea at that time (Berthon et al. 
2000 ). 

3.3.3 Size Distribution of Particles 

The size distribution of particles was measured using 
the Coulter Counter technique for the range 1.8-64 ^m. 
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Fig. 5. Time-series of a) total suspended matter concentration (Cj sm) 7 b) concentration of particles, No, 
from 2-30 ^m, and c) the coefficient £ of the Junge power law. The symbols and lines are as in Fig. 1. 


Figures 5b and 5c present the time series of the concentra- 
tion of particles, No (in units of per cubic centimeters) with 
diameters ranging from 2-30 /im and of the exponent £ of 
the Junge-type power function (Junge 1963) fitted using 
the least-squares method to the measured size distribution 
of particles: 

N = kD~S 9 (2) 

where N is the number of particles of diameter D, and k 
is the Junge coefficient. 

The concentration of particles, No (Fig. 5b), is well cor- 
related with CtsMi although a part of the measured Ctsm 
(0. 7-2.0 fim) is not included in the present size distribution 
of particles. Values of £ (Fig, 5c) range from 2. 5-5. 8 with 
an average of 4.06 (a = 0.62). Such a function does not al- 
ways adequately represent the measured size spectra, and 
probably causes some of the high variability in £ values 


observed for individual campaigns. Nevertheless, the use 
of this power function has ensured a first synthesis of the 
particle size distribution measurements. The derived high 
average value of £ is close to the one compiled by Jonasz 
(1983) for the Mediterranean Sea. The link between £ and 
the in situ attenuation and scattering coefficients will be 
discussed below. 

3.3.4 Vertical Distribution 

In the analyses presented above, as well as in those 
given later in this document, predominantly near-surface 
values were considered. The average seasonal profiles of 
temperature and salinity (Fig. 2) showed that the highest 
variability occurred within the near-surface layer (about 
0-5 m) in particular, because of the rapid passage of less 
salty water. Figure 6 shows the corresponding seasonal av- 
erage vertical profiles of C a and Ctsm* Two situations are 
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0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 2.5 

C a [mg rrT 3 ] CjSM^ m ” 3 ] 



Fig. 6. Average seasonal vertical profiles of C a , a) JFM, c) AMJ, e) JAS, and g) OND; and Ctsm, b) JFM, 
d) AMJ, f) JAS, and h) OND. The values are computed from the samples taken at three discrete depths 
(near-surface, 8m, and 14m). Standard deviations around the average are represented by dotted lines. 
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observed: 1) a rather homogeneous profile (e.g., spring or 
summer for C a ); and 2) stratification at the surface (e.g., 
winter for C a and Ctsm)- As seen with temperature and 
salinity, the near-surface layer shows the greatest variabil- 
ity Note that the shape of the Ctsm profile (and to a lesser 
extent the C a profile) is inversely correlated with respect 
to the average salinity vertical profile. 

The present vertical distributions of C a are in rather 
good agreement with the average seasonal distributions 
presented by Zavaterelli et al. (1998) for the shallow north- 
ern Adriatic Sea (north of the Po River estuary), although 
their values for the winter season were somewhat higher 
(average near-surface values of 2.7 mg m -3 ). In addition, 
the vertical profiles of the beam attenuation and absorp- 
tion coefficients (obtained from AC-9 measurements) de- 
monstrated that in almost all situations, the sampling at 
the three fixed depths (near-surface, 8 m, and 14 m) did 
not miss any particular aspect of the subsurface structure. 

3.4 Marine Inherent Optical Properties 

This section describes the variations of the measured 
inherent optical properties: in vivo absorption coefficients 
measured on filtered seawater samples, in situ absorption, 
and scattering coefficients estimated from the AC-9 mea- 
surements. 

3.4.1 In Vivo Absorption 

Spectra of the in vivo absorption coefficient of total and 
nonpigmented particulate matter, a p (A) and a<* p ( A), re- 
spectively, were measured within the range of 400-750 nm 
with a 1 nm resolution according to Tassan and Ferrari 
(1995) on the water samples collected at discrete depths. 
The difference of a p (A) — a<f p (A) gives the absorption co- 
efficient of pigmented matter, a p / i (A). The absorption co- 
efficient spectrum of colored dissolved organic matter, i.e., 
yellow substance, a ys ( A), was measured within the interval 
350-750 nm according to Ferrari et al. (1996). The spec- 
tra a p / l (A), ad p (A), and a y5 (A) are presented for the near- 
surface samples in Fig. 7 through their values at 443, 400, 
and 400 nm, respectively (statistics are given in Table 3). 
The values presented for ad p ( A) and a y5 (A) at 400 nm actu- 
ally are the ao values resulting from the following exponen- 
tial law fitted to each spectra (for the range 400-650 nm 
for a^p and 350-600 nm for a ys ): 

^xx(A) — noxx exp[ S XI (A Ao)] H- cb X xi (3) 

where xx is either dp or ys ; ao X x> S XX} and cb xx are the 
coefficient, slope, and background of the exponential law, 
respectively; and Aq is the reference wavelength (i.e., 400 
and 350 nm for dp and ys , respectively). 

The presence of cb xx in the fitting function deserves 
additional comment. Concerning a ya (A), it is well known 
that absorption in the red end of the spectrum is negligible 


(Kirk 1994). Measured absorption in this spectral range 
has been interpreted as a residual due to scattering by very 
small particles which passed through the (0.22 pm) filter, 
formed in the cuvette before the measurements (e.g., col- 
loids), or as a result of a difference in the refractive index 
between the sample and the reference due to a difference in 
salinity (Green and Blough 1994). In the results presented 
here, however, the background cb y3 exhibited an apprecia- 
ble variability for sampled water with the same salinity, 
therefore not supporting the Green and Blough (1994) hy- 
pothesis. On the contrary, a rough correlation between 
cb ys and a y5 (A) was observed, which supports the scatter- 
ing hypothesis. A constant background implicitly assumes 
scattering by submicron particles is not wavelength depen- 
dent, but this was chosen as such for the sake of simplicity. 

The same hypothesis of a null a^ p (750) value is not as 
well grounded, especially for inorganic particles. Never- 
theless, the general approach (Mitchell and Kiefer 1988) 
consists of subtracting the value measured at 750 nm from 
the whole spectra before fitting it to the exponential law. 
After several trials, (3) was found to give the best fit for the 
&dp spectra. In particular, it helped account for a change 
of slope frequently observed in the spectral domain 550- 
600 nm. Statistics on S xx and cb xx are given in Table 3. 

The a p h( 443) values range from 0.009-0.295 m -1 with 
an average of 0.0569 m -1 (a = 0.0391). As expected, the 
variations of ^^(443) closely followed the evolution of C a 
at seasonal (with spring and autumn maxima, winter and 
summer minima) and campaign time scales (Fig. 7a). The 
C a specific-absorption coefficient of the pigmented particu- 
late matter, a* h at 443 nm varies from 0.018-0.141 m 2 mg -1 
of C a , i.e., by about a factor of 8, with an average of 
0.0506 (a~ 0.0175). Its seasonal variation (Fig. 7b) is in- 
versely distributed compared to a p ^(443) because 
decreases when C a increases (e.g., Bricaud et al. 1995). 
Fig. 7e also presents the ratio a p / l (A):a p / l (443) for differ- 
ent SeaWiFS wavelengths. Variations of the presented ra- 
tios are rather limited. Roughly, low values of both ra- 
tios a p h (550)^^(443) and a p h(665):a p h(443) are associ- 
ated with high values of a p / l (488):a p /j(443) and vice-versa. 
Minimum values for the two first ratios may indicate a 
flattening of the absorption spectra due to the package 
effect (Morel and Bricaud 1981). In effect, minimal val- 
ues of a* h ( 443) are generally (but not systematically) ob- 
served in conjunction with these low ratios. The ratio 
a p / l (412):a p / l (443) shows some peaks with values higher 
than 1 probably due to high concentrations of phaeopig- 
ments (however this has not been confirmed by the HPLC 
measurements). 

The fitted a^ p (400) values range from 0.005-0.312 m" 1 
with an average of 0.0481 m -1 (a = 0.0318), whereas the 
average background cbd p value is 0.007 m -1 , i.e., about 
15% of the average a^ p (400) value. The seasonal evolu- 
tion of a<f p (400) (Fig. 7c) roughly follows that of a p h{ 443) 
(the plot of the one versus the other — not shown here — 
nevertheless exhibits a considerable scattering) and fol- 
lows quite well the evolution of the Ctsm concentration 
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Fig. 7. Time series, for near-surface samples, of a) the absorption coefficient of pigmented particulate matter 
a p h, at 443 nm in per unit meters; b) the corresponding chlorophyll Orspecific absorption coefficient a*/, (443) 
in meters square per milligram of chlorophyll a; c) the absorption coefficient of nonpigmented particulate 
matter aj p at 400 nm in per unit meters; and d) the absorption coefficient of colored dissolved organic matter 
a ys , at 400 nm in per unit meters. The symbols and lines are as in Fig. 1. 
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OCT95 JAN96 APR JUL OCT JAN97 APR JUL OCT JAN98 APR JUL OCTDEC98 


Fig. 7. (cont.) Panel e) shows the time series of the ratio a p /,(A)/a p / I (443) for A = 412 (solid line), 488 
(dotted line), 510 (dashed line), and 665 (dotted and dashed line). The values correspond to the averages of 
each campaign. 


Table 3. Surface absorption coefficients, in units of per meter, of pigmented ( a ph ) and nonpigmented (a dp ) 
particulate matter and of colored dissolved organic matter ( a ys ) at 443, 400, and 400 nm, respectively. Parameter 
(S xx and cb xx ) of the exponential curve fitted to the a dp and a ys spectra (see text). The period and statistic 
notations are the same as in Table 1. 


Period 

Statistic 

a p h(443) 

aj p (400) 

Sdp 


a ys (400) 

Sys 

c5 y ^ 

Global 


268 

268 

268 

268 

268 

268 

268 


X 

0.05685 

0.04807 

■ 



0.0181 

0.02605 


X„ 

0.03905 

0.03177 

| 

1 


0.0023 

0.03441 


X- 

0.0085 

0.0054 

■ 

1 

1 


-0.0343 


X* 

0.2954 

0.3119 





0.1969 

JFM 

X n 

64 

64 

64 

64 

64 

64 

64 


X 

0.06340 

0.04734 

0.0113 

0.00754 

0.11491 

0.0170 

0.02009 


X a 

0.06255 

0.02768 

0.0014 

0.00401 

0.03797 

0.0019 

0.03200 


X- 

0.0085 

0.0123 

0.008 



KEEK 

-0.0125 


A + 

0.2954 

0.1362 

0.015 


IB 

■ 

0.1603 

AMJ 

A n 

68 

68 

68 

68 

67 

67 

67 


X 

0.05601 

0.04579 

0.0126 

0.00746 

0.16232 

0.0190 

0.02794 


X a 

0.02987 

0.02960 

0.0013 

0.00556 

0.05380 

0.0021 

0.03523 


x~ 

0.0147 

0.0054 

0.010 

0.0000 

0.0666 

0.013 

-0.0343 


X* 

0.1671 

0.1602 

0.017 

0.0297 

0.3180 

0.026 

0.1969 

JAS 

X n 

62 

62 

62 

62 

63 

63 

63 


X 

0.04545 

0.04216 

0.0121 

0.00618 

0.14741 

0.0190 

0.03002 


Xa 

0.02428 

0.04183 

0.0010 

0.00574 

0.06186 

0.0029 

0.03326 


X- 

0.0170 

0.0090 

0.009 

-0.0005 

0.0482 

0.013 

-0.0116 


x + 

0.1549 

0.3119 

0.015 

0.0375 

0.3733 

0.028 

0.1643 

OND 

x n 

74 

74 

74 

74 

74 

74 

74 


X 

0.06153 

0.05576 

0.0115 

0.00709 

0.16798 

0.0176 

0.02613 


Xa 

0.02683 

0.02588 

0.0013 

0.00285 

0.05358 

0.0016 

0.03656 


X- 

0.0105 

0.0098 

0.008 

0.0021 

0.0791 

0.013 

-0.0256 


x + 

0.1325 

0.1117 

0.014 

0.0144 

0.2964 

0.022 

0.1664 


in Fig. 5a. The slope Sdp varies from 0.008-0.017 with an 
average of 0.0119 (cr = 0.0014), and the amplitude of the 
seasonal variations are lower than for ad p (Table 3). The 
slope Sdp does not exhibit any particular correlation with 

& dp * 

The a ya (400) values vary from 0.05-0.30 m _1 with a 
global average of 0.142 m -1 (a = 0.051) (Table 3). Here 
too, a global cycle with winter and summer minima and 


spring and autumn maxima can be identified although it is 
less evident for the period of April 1996-September 1996 
(Fig. 7d). The slope S ys in (3), resulting from the fit within 
the interval 350-600 nm, shows an average value of 0.0180 
(a = 0.0024). The resulting average background value cb ys 
is about 17% of the average a y5 (400) value. This back- 
ground may occasionally show negative values (minimum 
values in Table 3) that remain within the noise level. 


14 








J-F. Berthon, G. Zibordi, J.P. Doyle, S. Grossi, D. van der Linde, and C. Targa 



39 



Fig. 8. Time-series of a) the near-surface total (except water) scattering coefficient b t ~ w { A) at 555 nm 
(average value for aim layer); and b) the coefficient nb of the power function describing frt-u;(A) spectral 
dependency (see text). The symbols and lines are as in Fig. 1. 


3.4.2 In Situ Absorption and Scattering 

In situ vertical profiles of total (except water) spec- 
tral beam attenuation and absorption coefficients, c t — w (A) 
and a t _u,(A), respectively, at the nominal wavelengths 412, 
440, 488, 510, 555, 630, 650, 676, and 715 nm, were mea- 
sured with an AC-9 (Wetlabs, Inc.). The difference of 
c t _ w (A) — at-^A) gives the (non- water) spectral scattering 
coefficient b t - w ( A). The AC“9 deployments only started 
during campaign 21 (January 1997). 

In Zibordi et al. (2002), it was shown that the agree- 
ment between a t _ tlJ (A), corrected for scattering using Zan- 
eveld et al. (1992), and the sum of the in vivo absorp- 
tion coefficients of the three components described above 
[<Vi(A)+ad p (A)-ba y5 (A)] is reasonably good. This is partic- 
ularly true for the blue part of the spectrum, although the 
spreading of the points is rather high. Figure 8 presents the 
time series of b t - w {X) at 555 nm, i.e., c t _ lt ,(555)-a(_ u; (555) 
near the surface (average value within the 0.5-1. 5 m near- 
surface layer), together with the exponent nb of the power 
function describing its spectral dependence. This expo- 
nent was obtained by fitting each &t_tu(A) spectra within 
the interval 412-650 nm, to the following relationship 

b t - w (\) = b t _ w ( A 0 )[A/A 0 ] n6 , (4) 

where Ao = 555 nm. Quality control of the fit was per- 
formed by rejecting situations with a value exceeding a 
fixed threshold (1.10“ 3 ). 

Values of b t _ w ( 555) vary from 0.1-3. 5m -1 (Fig. 8a) 
with an average value of 0.96 m -1 (a — 0.68). Variations 


of 6 £ _u;( 555) are related to those of Ctsm (Fig. 5a) and of 
the concentration of particles N 0 (Fig. 5b). The values of 
the exponent nb shown in Fig. 8b vary from -0.9 to +0.2 
(actually, very few values are positive) with an average of 
—0.37 (a = 0.21). An approximate cycle can be identified 
with high values generally observed in the winter-spring 
period and lower ones in summer (this was true primar- 
ily for 1997). It is generally accepted that the spectra of 
the scattering coefficient become less flat as the particle 
populations are dominated by small particles (Stramski 
and Kiefer 1991). In winter-spring, the living particles are 
probably dominated by diatoms (this also corresponds to 
the period when river discharge may have a major influ- 
ence); whereas in summer, small species (e.g., cyanobacte- 
ria) are predominant, as observed in the pigment succes- 
sion presented in Sect. 3.3.1. 

Suspensions having a size distribution described by the 
Junge power law show a spectral dependence of the beam 
attenuation coefficient c( A) proportional to A 3 ~*, with £ 
being the exponent of the Junge law (Diehl and Haardt 
1980, and Reuter 1980). In the present situation, the aver- 
age value of 4.06 for £ would yield an average dependence 
of c t _u;(A) proportional to A 3-4 06 = A -106 . When di- 
rectly computed on measured individual Ct_ru(A) spectra 
within the interval of 412-650 nm, an average dependence 
with A-°- 81 ( ±0 - 23 ) was found. The value of the exponent 
lays between the typical —1 assumed for oceanic waters 
by Loisel and Morel (1998) and the mean -0.5 proposed 
by Voss (1992) for a wide range of oceanic and coastal 
locations. It also has to be recalled that the measured 
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size distributions were limited to 2 jim at the lower end of 
the diameter range and that many of them were not fully 
following the Junge law. According to Diehl and Haardt 
(1980) the expression c(A) oc A 3- ^ can also be applied to 
the wavelength dependence of the scattering coefficient. 
Nevertheless, as shown above, the spectral dependence of 
b t - w ( A) computed here is A -0 37 ( ±0 ' 21 \ suggesting an ap- 
preciable effect ofa t _ w (A) (particularly in the blue part of 
the spectrum) in the resulting average spectral dependence 
of ct_ w (A). 

Linear relationships relating b t - w ( A) and ^.^(555), as 
presented by Gould et al. (1999), were also investigated. 
The resulting slopes were close to those proposed by these 
authors, e.g., 1.063 at 412 nm and 0.869 at 676 nm versus 
1.160 and 0.850, respectively, given in Gould et al. (1999). 
Note that in this last analysis, the scattering contribution 
by pure water was included in the total scattering coeffi- 
cient, 6(A). 

3.5 Marine Apparent Optical Properties 

The marine apparent optical properties are derived from 
measured profiles of downward spectral irradiance Ed(z , A), 
upwelling spectral radiance L u (z, A), and upward irradi- 
ance E u (z t A). All the data are collected at the nominal 
center wavelengths of 412, 443, 490, 510, 555, 665, and 
683 nm using Satlantic ocean color irradiance- and radi- 
ance series-200 (OCT200 and OCR-200, respectively) ra- 
diometers. There are four variables derived from these 
measurements that are of particular interest: 

1) The diffuse attenuation coefficient for downward ir- 
radiance, K Ed { A), computed for an optically ho- 
mogenous near-surface layer (Zibordi et al. 2002), 
for which the minimum and maximum depths were 
typically 0.3 and 4.5 m, respectively; 

2) The subsurface (z = 0") remote sensing reflectance 
J? rs ( O’, A) = Z/ u (0~, A)/£d(0“, A) and irradiance re- 
flectance i?(0",A) = E u (0~ , X) / Ed(0~ y X); 

3) The factor of non-isotropy of the light field mea- 
sured at nadir, Qn(0“,A) = i?u(0“, A)/L u (0~, A); 
and 

4) The normalized water-leaving radiance, Lwn( A). 

As described in Zibordi et al. (2002), optical measure- 
ments and the resulting apparent optical properties were 
corrected for bottom effects (Maritorena et al. 1994), in- 
strument self-shading (Zibordi and Ferrari 1995) and tower 
shading (Doyle and Zibordi 2002). A quality control pro- 
cedure was performed on the data by eliminating the sta- 
tions for which any of the following were determined to be 
applicable: 

a. The ratio of diffuse to direct above-water down- 
ward irradiance was greater than 2 (in order to re- 


ject measurements not made under clear-sunf con- 
ditions); 

b. The relative difference of the irradiance reflectance 
(at 0“ depth) between two successive casts of the 
same station was greater than 10% (in order to re- 
ject measurements made in unstable} environmental 
conditions); or 

c. The K Ed ( A) spectral values computed for the near- 
surface layer were lower than the values measured 
by Smith and Baker (1978) for the clear waters of 
the Sargasso Sea. 

As a consequence, the number of optical stations used in 
the following analysis is reduced with respect to the prece- 
dent series. 

Figure 9 gives the time series of K Ed (5 55), Q„(0“, 555), 
and the ratio i?(0~, 490)/J?(0~, 555). From this point on, 
the depth argument z = 0’ is omitted for the sake of sim- 
plicity. The presentation and discussion on apparent prop- 
erties will be further developed in Sect. 4 which relates the 
seawater optical components and their associated inherent 
optical properties to the apparent optical properties. 

Figure 9a shows the diffuse attenuation coefficient for 
downward irradiance, K Ed at 555 nm, which varies from 
0.08-0.37m _1 , and has an average value of 0.180m" 1 
(<7=0.059). The tendency previously described for the opti- 
cal components (minimum in December- January and max- 
imum in March-April, as well as October-November) is 
roughly reproduced here for the years 1996 and 1997, how- 
ever, the year 1998 is much less defined. The K Ed { 555) 
values are well correlated with the Ctsm concentrations, 
but neither parameter exhibited a well-defined cycle during 
the year 1998 (Fig. 5a). 

The Q n values at 555 nm show amplitudes ranging from 
3.9-6.2sr with an average value of 4.87 (cr=0.49), which 
are in agreement with the theoretical studies of Morel and 
Gentili (1996). The temporal evolution of Q n show r s al- 
most regular large variations overlapped by relatively small 
short-term variations (Fig. 9b). The large variations, ex- 
hibiting a maximum in winter and a minimum in summer, 
are mainly attributed to seasonal changes of the average 
solar zenith angle. The short-term variations, within cam- 
paigns, are attributed to changes in the water optical prop- 
erties and in the illumination conditions (variations of the 
sun zenith, aerosol content, cloudiness, etc.) from station 
to station. 

The analyses of Q n spectra also show a significant de- 
pendence with wavelength. The results exhibit higher val- 
ues of Q n at wavelengths characterized by higher absorp- 
tion of seawater, mostly in the red part of the spectrum. In 


f Refer to Zibordi et al. (2002) for an in-depth discussion of 
the definition of clear-sun and clear-sky conditions, and the 
differences between the two. 

} Principally those data affected by illumination changes, sea 
state effects, etc., and as further explained in Zibordi et al. 
( 2002 ). 


16 



J-F. Berthon, G. Zibordi, J.P. Doyle, S. Grossi, D. van der Linde, and C. Targa 





OCT95 JAN96 APR JUL OCT JAN97 APR JUL OCT JAN98 APR JUL 0CTDEC98 


Fig. 9. The time-series’ of a) K Ed { 555) computed for an optically homogeneous near-surface layer; b) the 
surface (0") Q n factor at 555 nm, and c) the surface (O') reflectance ratio #(490)/i?(555). The symbols and 
lines are as in Fig. 1. 


the spectral range 412-665 nm, the highest average value 
(i.e., 5.56 with cr=0.99) and the highest annual variations 
(i.e., ranging from approximately 4.0-8. 2) have both been 
observed at 665 nm. A more detailed analysis on the Q- 
factor dependence with the sun geometry and the optical 
properties of the seawater is given in Zibordi and Berthon 
( 2001 ). 

Figure 9c presents the time series of the irradiance re- 
flectance ratio R(490)/R(555). The two wavelengths for 
this ratio are classically coupled when trying to empirically 
relate reflectance or normalized water-leaving radiance ra- 
tios to the chlorophyll a concentration in the open ocean 
(e.g., O’Reilly et al. 1998 and Kahru and Mitchell 1999). 
Values of this ratio encompass a rather limited range, from 
0. 5-2.0, with an average value of 1.08 ( <7=0. 27). The vari- 
ability within individual campaigns is small compared to 
the previously studied quantities. It can clearly be seen 


that R(490)/R(555) is inversely correlated with i^£ d (555) 
shown in Fig. 9a. The presentation of relationships be- 
tween reflectance ratios and optical components will be 
made in the following sections. 

4. DISCUSSION 

In this section, the bio-optical status of the site (i.e., 
its classification according to Case-1 or Case-2 types) is 
discussed and empirical algorithms for chlorophyll a and 
total suspended matter retrieval from marine reflectance, 
are presented. 

4.1 Case-1 and Case-2 Waters 

The Case-1 water corresponds to water for which the 
optical properties are only determined by the phytoplank- 
ton and its derivative products (Morel and Prieur 1977). 
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On the contrary, the optica] properties of the Case-2 waters 
are also determined by other material (from terrestrial- or 
bottom origin for example) not systematically co-varying 
with the phytoplanktonic components. 

In the present study, the discrimination between Case- 
1 and Case-2 water types is based on the recent work of 
Loisel and Morel (1998), who revisited the relationship 
between C a and the total (except water) beam attenu- 
ation coefficient at 660 nm, c t - w (660). Here, the values 
of c t - w {660) have not been directly measured — they have 
been extrapolated from the measurements at 670 nm us- 
ing the average slope of the Cf_^(A) spectral dependence 
presented above. The average value computed within the 
0.5-1. 5 m near-surface layer is related to the near-surface 
C a value (Fig. 10). According to the discrimination line 
proposed by Loisel and Morel (1998), about 63% of the 
CoASTS stations revealed a Case-1 behavior (dashed line). 
Note that this differentiation has been possible for only 
about half of the data set, that is, when in situ measure- 
ments of c f _ ti; (A) were available. 



C G [mg m -3 ] 


Fig. 10. The relationship, for the near-surface 
layer, between C a and the total (except water) beam 
attenuation coefficient c^^A) at 660 nm (average 
value for aim layer). The solid line represents the 
linear regression computed for the present data set. 
The dashed line represents the limit proposed by 
Loisel and Morel (1998) segregating Case-1 (open 
circles) and Case-2 (solid circles). The dotted line 
is the limit adopted by Sagan (1995). 

Figure 10 also presents a relationship, c t - w (660) = 
0.5C a , of the same type adopted by Mitchell and Holm- 
Hansen (1991) and Sagan et al. (1995) for classifying Case- 
1 and Case-2 waters (dotted line) in the Antarctic Penin- 
sula area. Using this latter threshold, the percentage of 
Case-1 stations for the present study was only 40.5%. The 
divergence between the two relationships is of particular 
importance within the C a interval 0. 1-2.0 mg m~ 3 . This 


illustrates the difficulty in segregating Case-1 and Case-2 
waters without ambiguity on the basis of a simple thresh- 
old in coastal areas such as the platform site. 

4.2 Empirical Algorithms 

The potential of this data set for ocean color remote 
sensing applications in the northern Adriatic Sea was eval- 
uated by examining empirical relationships between sur- 
face reflectance and the chlorophyll a concentration. Rela- 
tionships envisaged here involve the in-water remote sens- 
ing reflectance, R TS ( A) = L u (\)/E d (\). As an example, 
Fig. 11 presents the global relationship for the CoASTS 
data between C a (in milligrams per cubic meter) and the 
reflectance ratio R rs (490)/R rs (555), with the latter denoted 
in the more compact form for future formulations. As 
a way of comparison with the SeaWiFS Ocean Color 2 ver- 
sion 4 (OC2v4) algorithm (O’Reilly et al. 2000), a third- 
order polynomial was fit to the log-transformed (base 10) 
data and the resulting relationship is: 

log 10 [a] = 0.091 - 2.620 log 10 [i£ 5 ] 

- 1.148 1og 10 [/Z r 3 s 5 ] 2 (6) 

- 4.949 log 10 [i? 3 ' 5 ] 3 . 
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Fig. 11 . Relationships between the remote sensing 
reflectance ratio (0‘ depth) R rs (490)/ J R rs (555) and 
C a . Filled and empty circles represent the Case- 
2 and Case-1 situations, respectively (according to 
Loisel and Morel 1998), whereas the stars repre- 
sent the unidentified situations. The thick solid line 
represents the polynomial fit (see text), the thick 
dotted-dashed line the linear fit. The dotted and 
dashed thin lines represent a factor of ±1.5 and ±2 
around the fit, respectively. The thick dashed line 
is the SeaWiFS OC2v4 algorithm. 
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The global relationship for the CoASTS data set given 
in (6) explains approximately 77% of the variance in C a 
and the root mean square of relative difference (RMSrd) 
is approximately 42%. In 34 situations out of 109 (i.e., in 
about 31% of the cases), the fitted C a values differ from 
the measured values by more than a factor of 1.5. The 
OC2v4 algorithm (thick dashed line in Fig. 11) systemat- 
ically overestimates C a when applied to the present data 
set, particularly in the range of 0. 1-1.0 mg m” 3 . For equal 
C a values, the corresponding reflectance ratio is lower than 
that given by the OC2v4 algorithm. The overestimate for 
the latter is probably because of the increased absorption 
by dissolved and particulate matter in the blue part of 
the spectrum. The agreement is much better in the range 
of 2.0-10.0 mg m -3 when C a is probably the dominating 
optical component. 

For the purposes of comparison with (6), a linear re- 
lationship between log-transformed chlorophyll a and re- 
flectance values was computed and gives very similar per- 
formances (where the determination coefficient r 2 equals 
0.76 and RMSrd = 0.43): 



^TSM [9 m 1 


Fig. 12. Relationships between the remote sens- 
ing reflectance ratio (0" depth) i? rB (510)/R ra (665) 
and Ctsm- The solid line represents the linear fit 
(see text). The thin dotted and dashed lines and 
symbols are as in Fig. 11. 


log 10 [C.] = 0.079 - 2.898 log 10 [it* 5 ]. (7) 

The two relationships behave very closely within the range 
0.3-10 mg m" 3 and start to diverge for C a values lower 
than 0.3 mg m -3 (Fig. 11). The application of these algo- 
rithms outside the range used for their development, how- 
ever, in particular for C a values lower than O.lmgm* 3 , 
should be considered with caution. 

Note that the use of the subsurface remote sensing 
reflectance instead of the above-water remote sensing re- 
flectance in such algorithms, does not make an appreciable 
difference when reflectances are used in ratios. Among the 
subset of stations for which it was possible (i.e., when AC-9 
data were available together with C a values), the applica- 
tion of the Loisel and Morel (1998) relationship resulted 
in 49 and 51% of Case-1 (open circles) and Case-2 (solid 
circles) situations (Fig. 11), respectively. 

Figure 12 presents the global relationship between val- 
ues of Ctsm (in grams per cubic meter) and the in-water 
remote sensing reflectance ratio, R rs (510)/i? rs (665), corre- 
sponding to the following log-linear fit: 


iogioI^TSM] = 0.967 


1.169 


log 10 


flrs(510) 

R r8 (665)_ * 


( 8 ) 


This relationship explains 66% of the variance in Ctsm, 
with a RMSrd of 39%. In 31 situations out of 109 (about 
28 %), the fitted Ctsm values differ from the measured val- 
ues by more than a factor of 1.5. 


The use of other relationships involving more than two 
reflectances at two different channels was investigated. In 
particular, multilinear regressions between C a or Ctsm 
and several reflectances (up to 6) sometimes explained a 
higher percentage of variance. These relationships, how- 
ever, have been rejected because of the evidence of high 
intercorrelation among reflectances, which can lead to a 
decrease in the predictive power of the relationships. Fig- 
ure 13 presents the relationship between the ratio of nor- 
malized water-leaving radiance Lvw(490)/Lwtv(555) and 
the diffuse attenuation coefficient for downward irradiance, 
^£*(490), computed for an optically homogeneous near- 
surface layer. The corresponding local algorithm, equiv- 
alent to the SeaWiFS algorithm (Mueller 2000) for open 
ocean application, is the following: 


#^(490) = 0.016 + 0.205 


L\vn{ 490) 

Lwat(555) 


- 1.754 


( 9 ) 


where 0.016 m -1 is the constant value of K Ed (490) for pure 
water (adopted by Mueller 2000). This relationship ex- 
plains 80% of the variance in K Ed { 490), with an RMSrd of 
23%. 

The relatively high scattering of the points in Fig. 13 
within the K Ed interval 0.1-0.5 m* -1 (leading to an r 2 value 
of only 0.80) may be partly explained by: 

■ The variability in the relationships between C a and 
Ctsm (yellow substance influence is relatively lower 
at these wavelengths) , in particular for high concen- 
trations, inducing changes in the radiance spectral 
shape; and 
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■ The uncertainty in the retrieval of Lwn and K£ d 
in situations of strong vertical (subsurface) stratifi- 
cation. 

The slope of the present algorithm is slightly higher than 
for the SeaWiFS algorithm (—1.5401) proposed by Mueller 
(2000) and as such, leads to significantly higher values of 
KE d for the range 0.2-0.5m _1 . [Actually, the maximum 
value of KE d { 490) was 0.25 m -1 within the data set used 
to build the SeaWiFS algorithm.] 



0.01 0.10 1.00 
K Ed (490) [m-'] 


Fig. 13. The relationship between the normalized 
water-leaving radiance ratio Lvvn(490)/Lvvn(555) 
and the diffuse attenuation coefficient Ks d { 490). 
The solid line represents the fitted power function, 
the dotted and dashed lines represents a factor of 
±1.5 and ±2 around the fit, respectively. The dot- 
ted and dashed line represents the Mueller (2000) 
algorithm. The symbols are as in Fig. 11. 

5. CONCLUSIONS 

The first three years (from autumn 1995 to autumn 
1998) of the CoASTS data collected from an oceanographic 
tower in the northern Adriatic Sea have been presented 
and analyzed. The high degree of completeness of the 
data set, including both marine and atmospheric measure- 
ments, makes it appropriate for calibration and validation 
activities associated with ocean color sensors (Sturm and 
Zibordi 2002). Care must be taken with the inevitable 
time lag when creating the matchups between the satellite 
imagery and the in situ measurements, because variability 
on the scale of one half to one full day cannot be neglected 
at this site. For the marine components, this short-term 
variability — mainly affecting the near-surface layer — is su- 
perimposed on a classical seasonal variability and is well 
correlated with salinity variations at the surface resulting 
from the southward transport of water discharged by rivers 
into the northern Adriatic Sea. 


The concentrations of chlorophyll a ( C a ) and total sus- 
pended matter (Ctsm) span over two orders of magnitude. 
A similar range of expression is observed for the absorption 
coefficients of particulate and dissolved material. Absorp- 
tion by pigmented and nonpigmented particulate matter 
are well correlated with C a and Ctsm> respectively. In 
addition, the spectral dependence of the absorption coef- 
ficients exhibits an appreciable variability (e.g., a factor 2 
for the slopes describing the exponential decrease from the 
blue to the red spectral domain of nonpigmented particu- 
late and of dissolved matter absorption coefficients). Al- 
though no measurement of the back-scattering coefficient 
was undertaken up to the year 2000, the measured particle 
size distribution could be used within Mie computations 
to predict the scattering phase function of particles and, 
hence, the backscattering coefficient. The coefficient of 
the Junge power law fitted to these size distributions, is in 
agreement with values previously found for the Mediter- 
ranean Sea and also shows a seasonal variability. 

The site presents bio-optical characteristics pertaining 
to both Case-1 (about 63% of the time) and Case-2 (ap- 
proximately 37% of the time) waters according to the clas- 
sification relationships between the beam attenuation coef- 
ficient at 660 nm, c t - w ( 660), and the C u proposed by Loisel 
and Morel (1998). Nevertheless, it was shown that the clas- 
sification of coastal areas, such as the platform site, based 
on a single unique segregation threshold is very challeng- 
ing. The measured apparent properties (Aj 5 d (A), Q n (A), 
R(A), etc.] also encompass a rather wide range of values; 
their variability is in agreement with the cycle described 
for the optical components and their inherent optical prop- 
erties. It has been highlighted that the Q-factor (measured 
at nadir) shows a significant spectral dependence and that 
short-term variations, in particular those due to changes 
in the water content, are superimposed over the seasonal 
cycle (mostly defined by the average sun zenith cycle). 

The local empirical algorithm making use of the ratio 
i? rs (490)/ii r8 (555) explains 77% of the variance in the C a 
near-surface concentration and is likely to provide lower 
values (by a factor of 2-4) in the range of 0. 1-1.0 mg m -3 
than the SeaWiFS OC2v4 algorithm, whereas the two algo- 
rithms converge in the range 1.0-10.0 mg m~ 3 . A local em- 
pirical algorithm making use of the ratio J R rs (510)/H rs (665) 
explains 66% of the variance in the Ctsm near-surface con- 
centration. The algorithm using the ratio of normalized 
water-leaving radiances Twn(490)/Lvvw(555) explains 80% 
of the variance in ifg d (490). The rather large range of C 0 
values covered by the data set probably makes the C a algo- 
rithm applicable to the whole northern Adriatic Sea with 
the possible exception of some unique situations, such as, 
the area very close to the coast or in near proximity to 
river plumes (e.g., the Po River plume). 
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Glossary 

AAOT Acqua Alta Oceanographic Tower 
AERONET Aerosol Robotic Network 
AMJ April-May-June 
AMT Atlantic Meridional Transect 

CoASTS Coastal Atmosphere and Sea Time Series 
CTD Conductivity, Temperature, and Depth 
CZCS Coastal Zone Color Scanner 
GSFC Goddard Space Flight Center 
HPLC High Performance Liquid Chromatography 

IOCCG International Ocean Colour Coordinating 
Group 

JGOFS Joint Global Ocean Flux Study 
JAS July-August-September 
JFM January-February-March 
MERIS Medium Resolution Imaging Spectrometer 
MOBY Marine Optical Buoy 

MODIS Moderate Resolution Imaging Spectroradiome- 
ter 

MOS Modular Opto-electronic Sensor 

NAd North Adriatic (Current) 

NASA National Aeronautics and Space Administra- 
tion 

OC2v4 Ocean Chlorophyll 2 (algorithm) version 4 
OCI Ocean Color Irradiance 
OCR Ocean Color Radiance 
OCTS Ocean Color and Temperature Scanner 
OND October-November-December 

PlyMBODy Plymouth Marine Bio-Optical Data Buoy 
PSU Practical Salinity Units 

RMSrd Root Mean Square of relative difference 
SeaWiFS Sea- viewing Wide Field-of-view Sensor 
TSM Total Suspended Matter 

UNESCO United Nations Educational, Scientific, and 
Cultural Organization 

YBOM Yamato Bank Optical Mooring 
Symbols 

Udp(A) Nonpigmented particulate matter absorption coeffi- 
cient. 

a p (A) Total particulate matter absorption coefficient. 
a p /v(A) Pigmented particulate matter absorption coefficient. 
a p/i(A) Chorophyll Or-specific pigmented particulate matter 
absorption coefficient. 


at-™ (A) Total (except water) absorption coefficient. 
a y5 (A) Colored dissolved organic matter (yellow substance) 
absorption coefficient. 
a yxx Generic symbol for aod P or ao yJ . 
ao dp Coefficient of a,dp( A) exponential law. 
aoys Coefficient of a ya (A) exponential law. 

6(A) Total scattering coefficient. 

6 t -™(A) Total (except water) scattering coefficient. 

c(A) Total beam attenuation coefficient. 

Ct-™(A) Total (except water) beam attenuation coefficient. 
cbdp Background of ad p (A) exponential law. 
cby 3 Background of a ya (A) exponential law. 
cb xx Either cbdp or cb y3 . 

C Pigment concentration. 

C a Chlorophyll a concentration. 

Ca Alloxanthin concentration. 

Cb Chlorophyll 6 concentration. 

Cd Diatoxanthin concentration. 

Cp Diadinoxanthin concentration. 

Cf Fucoxanthin concentration. 

Chb The sum of the 19'-hexanoyloxy fucoxanthin and 19'- 
butanoyloxy fucoxanthin concentrations. 

Ctsm Total suspended matter concentration. 

Cz Zeaxanthin concentration. 

Ci 2 Chlorophyll ci + C 2 concentration. 

D Particle diameter. 

Ed(z } X) Downward irradiance at depth z. 

Ed( 0”,A) Subsurface downward irradiance. 

E u (z y A) Upward irradiance at depth z. 

E u ( CT,A) Subsurface upward irradiance. 
k Junge coefficient. 

KE d ( A) Diffuse attenuation coefficient of Ed{z, A). 

L u (zy A) Upwelling radiance at depth z. 

L u ( 0", A) Subsurface upwelling radiance. 

Lvwv(A) Normalized water leaving radiance, 
m Air mass, 
n Number of samples. 

N Number of particles of diameter D. 

No Concentration of particles. 
nb Exponent of the b t - w {X) power function. 

P a Atmospheric pressure. 

<3n(0~,A) Q - factor for nadir, £u(0 , X)/Lu(0 ,A). 
r 2 Determination coefficient. 

R( 0”,A) In- water reflectance, E u (0 , A)/JEd(CT, A). 

R rs ( 0”,A) In- water remote sensing reflectance, defined as the 
L u (0~ , X) / E d {0~ } X) ratio. 

R r 3 B 5 The R rs (0",490)/ J R rs (0",555) ratio. 

RH Relative humidity. 

Sd P Slope of a<i p (A) exponential law. 

S w Seawater salinity. 

S xx Either Sdp oi* Sy 3 . 

S ys Slope of a ya (A) exponential law. 

T a Air temperature. 

T™ Water temperature. 

Wd Wind direction. 

W 3 Wind speed. 

xx Generic symbol for dp or ys. 

X The average value of variable X . 

X~ The minimum value of variable X. 

X + The maximum value of variable X. 

X n The number of samples for variable X. 

X 0 The standard deviation of variable X. 
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z Depth. 

a Angstrom exponent. 

0 Angstrom coefficient. 

7 Density excess. 

A Wavelength. 

Ao Reference wavelength. 

£ Junge exponent. 

< 7 Standard deviation. 
r 0 (A) Aerosol optical thickness. 

X 2 Chi squared. 
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remote sensing reflectance ratios and optical components or properties of interest (chlorophyll a y total suspended 
matter, and the diffuse attenuation coefficient) are presented. 
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